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ABSTRACT 

The  star  HD  193322  is  a  remarkable  multiple  system  of  massive  stars  that  lies  at  the  heart  of  the  cluster 
Collinder  419.  Here  we  report  on  new  spectroscopic  observations  and  radial  velocities  of  the  narrow-lined  component 
Abl  which  we  use  to  determine  its  orbital  motion  around  a  close  companion  Ab2  (P  —  312  days)  and  around  a 
distant  third  star  Aa  ( P  —  35  years).  We  have  also  obtained  long  baseline  interferometry  of  the  target  in  the  K' 
band  with  the  CHARA  Array  which  we  use  in  two  ways.  First,  we  combine  published  speckle  interferometric 
measurements  with  CHARA  separated  fringe  packet  measurements  to  improve  the  visual  orbit  for  the  wide  Aa,Ab 
binary.  Second,  we  use  measurements  of  the  fringe  packet  from  Aa  to  calibrate  the  visibility  of  the  fringes  of  the 
Abl,Ab2  binary,  and  we  analyze  these  fringe  visibilities  to  determine  the  visual  orbit  of  the  close  system.  The  two 
most  massive  stars,  Aa  and  Abl,  have  masses  of  approximately  21  and  23  M0,  respectively,  and  their  spectral  line 
broadening  indicates  that  they  represent  extremes  of  fast  and  slow  projected  rotational  velocity,  respectively. 
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1.  INTRODUCTION 

Massive  O-type  stars  are  usually  found  with  one  or  more 
nearby  companions  (Mason  et  al.  1998,  2009).  Most  of  these 
luminous  stars  are  very  distant,  and  consequently,  we  generally 
only  detect  their  very  nearby  companions  through  their  Doppler 
shifts  or  very  distant  companions  that  are  angularly  resolved. 
We  must  rely  on  high  angular  resolution  observations  of  the 
few  nearby  cases  to  detect  those  elusive,  mid-range  separation 
binary  stars.  One  of  the  most  revealing  examples  is  HD  193322 
(09  V:((n));  Walborn  1972),  the  central  star  in  the  sparse  open 
cluster  Collinder  419.  The  distance  to  the  cluster  is  741  ±  36  pc 
according  to  a  recent  study  by  Roberts  et  al.  (2010).  The  star’s 
complex  multiplicity  became  apparent  with  the  discovery  of 
a  companion  Ab  through  speckle  interferometry  observations 
by  McAlister  et  al.  (1987).  They  designated  the  system  as 
CHARA  96  Aa  (McAlister  et  al.  1989),  and  subsequent  speckle 
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measurements  detected  its  orbital  motion  (Hartkopf  et  al. 
1993;  Hartkopf  2010).  The  Aa,Ab  pair  was  also  recently 
resolved  through  the  technique  of  lucky  imaging  by  Malz 
Apellaniz  (2010).  The  composite  optical  spectrum  is  dominated 
by  a  relatively  narrow-lined  component  Abl,  and  Fullerton 
(1990)  discovered  significant  radial  variations  in  this  component 
indicative  of  a  spectroscopic  binary.  The  first  spectroscopic 
orbit  for  Abl  was  presented  by  McKibben  et  al.  (1998),  who 
determined  an  orbital  period  of  311  days.  In  addition  to  the 
close  Abl,Ab2  spectroscopic  pair  and  the  speckle  Aa,Ab  pair, 
there  is  another  wider  companion  B  at  an  angular  separation 
of  2.68  arcsec  (Turner  et  al.  2008).  Components  C  and  D  are 
more  distant  companions  that  also  occupy  the  central  region  of 
Collinder  419  (Roberts  et  al.  2010),  and  it  is  uncertain  whether 
they  are  orbitally  bound  to  the  central  multiple  system.  A  mobile 
diagram  presenting  the  known  components  of  the  system  is 
illustrated  in  Figure  1 .  The  long  orbital  period  estimate  for  A,B 
pair  is  based  upon  the  probable  masses  (see  Table  8  below), 
distance,  and  the  assumption  that  the  projected  separation  is  the 
semimajor  axis. 

Measurements  of  the  orbital  motions  of  the  stars  in  this  sys¬ 
tem  offer  us  the  means  to  estimate  the  masses  of  the  compo¬ 
nents.  We  have  continued  our  interferometric  and  spectroscopic 
monitoring  of  the  system  over  the  last  decade,  and  here  we 
present  a  progress  report  on  the  orbits,  mass  estimates,  and  spec¬ 
tral  properties  of  the  component  stars.  Combined  speckle  and 


1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

JUL  2011 


2.  REPORT  TYPE 


4.  TITLE  AND  SUBTITLE 

An  Interferometric  And  Spectroscopic  Analysis  Of  The  Multiple  Star 
System  HD  193322 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.  S.  Naval  Observatory, 3450  Massachusetts  Avenue, 
NW, Washington, DC, 20392 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED 

00-00-2011  to  00-00-2011 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROIECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  Astronomical  Journal,  142:21  (14pp),  2011  July 

14.  ABSTRACT 

The  star  HD  193322  is  a  remarkable  multiple  system  of  massive  stars  that  lies  at  the  heart  of  the  cluster 
Collinder  419.  Here  we  report  on  newspectroscopic  observations  and  radial  velocities  of  the  narrow-lined 
component  Abl  which  we  use  to  determine  its  orbital  motion  around  a  close  companion  Ab2  (P  =  312  days) 
and  around  a  distant  third  star  Aa  (P  =  35  years).  We  have  also  obtained  long  baseline  interferometry  of 
the  target  in  the  K  band  with  the  CHARA  Array  which  we  use  in  two  ways.  First,  we  combine  published 
speckle  interferometric  measurements  with  CHARA  separated  fringe  packet  measurements  to  improve  the 
visual  orbit  for  the  wide  Aa,Ab  binary.  Second,  we  use  measurements  of  the  fringe  packet  from  Aa  to 
calibrate  the  visibility  of  the  fringes  of  the  Abl,Ab2  binary,  and  we  analyze  these  fringe  visibilities  to 
determine  the  visual  orbit  of  the  close  system.  The  two  most  massive  stars,  Aa  and  Abl,  have  masses  of 
approximately  21  and  23M  ,  respectively,  and  their  spectral  line  broadening  indicates  that  they  represent 
extremes  of  fast  and  slow  projected  rotational  velocity,  respectively 

15.  SUBIECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

15 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


The  Astronomical  Journal,  142:21  (14pp),  201 1  July 


TEN  BRUMMELAAR  ET  AL. 


HD  193322 


Ab 


A 


BCD 


P  =  11000  y 
Aa 


P  =  35  y 


B1.5V  B8V  B9V 


Abl 


Ab2 

P  =  0.9  y 


09  Vnn 


08.5  111  [B2.5V] 

Figure  1.  Mobile  diagram  of  the  components  of  the  multiple  star  HD  193322. 
The  spectral  classification  is  given  under  each  stellar  component.  The  classifi¬ 
cation  for  Ab2  is  enclosed  in  brackets  to  emphasize  its  uncertainty  (based  upon 
its  relative  flux  contribution;  see  Sections  3  and  5).  The  classifications  for  C  and 
D  are  from  Roberts  et  al.  (2010).  The  period  estimate  for  A,B  is  based  upon  the 
projected  separation  (Mason  et  al.  1998). 


interferometric  observations  of  the  motion  of  the  Aa,Ab  pair  are 
used  in  Section  2  to  derive  a  preliminary  orbit  for  the  wide  sys¬ 
tem.  In  Section  3,  we  present  new  long  baseline  interferometric 
measurements  of  the  Abl,Ab2  pair  that  are  calibrated  using  the 
visibility  of  the  Aa  companion.  In  Section  4,  we  describe  a  di¬ 
verse  collection  of  spectroscopic  observations  which  we  use  to 
derive  a  revised  orbit  for  the  narrow-lined  Abl  component  in 
the  close  pair.  In  Section  5,  we  apply  a  Doppler  tomography 
algorithm  to  a  subset  of  the  blue  spectra  to  extract  the  spectra 
of  the  components.  Finally  in  Section  6,  we  discuss  the  masses 
and  other  properties  of  the  components  of  the  system. 

2.  VISUAL  ORBIT  OF  THE  WIDE  SYSTEM 

The  orbital  motion  of  the  Aa,Ab  pair  has  been  followed 
since  its  discovery  through  continued  speckle  interferometry 
observations  made  mainly  with  the  Mayall  4  m  telescope  at 
Kitt  Peak  National  Observatory  (KPNO;  McAlister  et  al.  1989, 
1993;  Mason  et  al.  1998,  2009).  The  date,  position  angle  9, 
and  separation  p  of  these  previously  published  observations  are 
collected  in  Table  1  for  convenience.  While  outliers  exist,  for  this 
magnitude  difference  Am  and  p  regime,  the  errors  from  speckle 
interferometry  measures  are  approximately  0.  5  in  position  angle 
and  0.5%  in  separation.  We  have  also  measured  the  relative 
motion  through  optical  long  baseline  interferometry  (OLBI) 
with  the  GSU  CHARA  Array  at  Mount  Wilson  Observatory 
(ten  Brummelaar  et  al.  2005).  These  separation  and  position 
angle  measurements  are  determined  by  measuring  the  fringe 
packet  separation  when  possible  along  two  pairs  of  baselines 
with  approximately  orthogonal  directions  projected  onto  the 
sky  (the  separated  fringe  packet  or  SFP  method;  Farrington  et  al. 
2010).  This  separation  is  determined  by  fringe  fitting  in  order 
to  avoid  shifts  caused  by  overlap.  Other  methods,  like  fits  of 
the  fringe  envelope,  may  suffer  if  one  fringe  packet  overlaps  the 
secondary  lobe  of  the  other  and  causes  the  center  of  the  fringe 
envelope  to  move.  Like  lunar  occultation  measurements,  a  single 
baseline  measurement  provides  a  separation  in  one  direction 
only.  Each  of  these  measurements  defines  a  line  in  astrometric 
space,  and  observations  at  several  projected  angles  are  required 
to  fully  define  the  position  of  the  secondary.  The  location  of  the 


Table  1 

Astrometric  Measurements  of  Aa,Ab 


Date 

e 

P 

Data 

Reference 

(BY) 

(deg) 

(arcsec) 

Type 

0.049  Speckle 

0.049  Speckle 


1985.5177 

188.4 

1985.8396 

192.5 

1986.8884 

198.6 

1988.6630 

216.6 

1989.7061 

229.6 

2005.6054 

109.1 

2005.7350 

107.7 

2005.8652“ 

100.4 

2006.4324“ 

100.1 

2006.4897 

101.5 

2006.5881 

113.9 

2006.6758 

118.0 

2007.4729 

111.7 

2007.5098 

113.6 

2007.6042“ 

100.9 

2008.4508 

116.8 

2008.6198 

121.8 

2008.8028 

124.7 

2009.4178 

120.1 

2009.5017 

126.7 

2009.6146 

122.0 

2009.7776 

122.0 

2010.8753 

129.8 

0.049 

Speckle 

0.048 

Speckle 

0.045 

Speckle 

0.0638 

OLBI/SFP 

0.0647 

OLBI/SFP 

0.086 

Speckle 

0.0409 

OLBI/SFP 

0.0670 

OLBI/SFP 

0.0651 

OLBI/SFP 

0.0565 

OLBI/SFP 

0.0666 

OLBI/SFP 

0.0665 

OLBI/SFP 

0.067 

Speckle 

0.066 

Speckle 

0.0616 

OLBI/SFP 

0.0551 

OLBI/SFP 

0.0626 

OLBI/SFP 

0.0575 

OLBI/SFP 

0.0651 

OLBI/SFP 

0.0649 

OLBI/SFP 

0.0648 

OLBI/SFP 

McAlister  et  al.  (1993) 
McAlister  et  al.  (1993) 
McAlister  et  al.  (1993) 
McAlister  et  al.  (1993) 
McAlister  et  al.  (1993) 
This  paper 
This  paper 
Mason  et  al.  (2009) 
This  paper 
This  paper 
This  paper 
This  paper 
This  paper 
This  paper 
Mason  et  al.  (2009) 
Mason  et  al.  (2009) 
This  paper 
This  paper 
This  paper 
This  paper 
This  paper 
This  paper 
This  paper 


Note.  a  Assigned  zero  weight  in  the  fit. 

(FITS  files  of  this  table  are  available  in  the  online  journal.) 


secondary  is  defined  as  the  point  with  the  minimum  total  rms 
distance  from  these  lines,  as  weighted  by  the  variance  of  the 
fringe  separations.  Formal  errors  are  calculated  using  a  method 
analogous  to  a  /  2  analysis,  and  the  errors  for  p  and  9  are  defined 
as  the  distance  change  required  to  increase  the  weighted  rms  by 
1.0  in  /2.  These  data  have  fairly  low  signal-to-noise  ratio  (S/N) 
and  for  many  epochs  we  have  data  from  only  a  single  baseline 
with  a  varying  position  angle  from  diurnal  motion,  so  there 
is  more  scatter  in  the  resulting  astrometry  than  in  the  speckle 
data.  All  the  speckle  and  CHARA  measurements  are  collected 
in  Table  1 . 

We  determined  a  new  orbital  solution  for  the  wide  pair  using 
the  combined  set  of  speckle  and  long  baseline  interferometric 
observations.  Note  that  we  did  not  include  the  measurement 
from  Mafz  Apellaniz  (2010)  using  lucky  imaging  with  the 
AstraLux  instrument  because  of  its  relatively  large  error.  All 
the  measurements  were  initially  assigned  equal  weight,  but  in 
the  orbital  fitting  process  we  identified  three  discrepant  points 
with  large  residuals  that  we  subsequently  zero  weighted  in  the 
final  fit  (those  dates  are  marked  in  Table  1).  The  orbit  was 
determined  using  the  grid  search  method  described  by  Hartkopf 
et  al.  (1989).  The  orbital  elements  are  listed  in  Table  2  and  the 
appearance  of  the  visual  orbit  on  the  sky  is  shown  in  Figure  2. 
The  original  speckle  data  set  covers  about  12%  and  the  recent 
speckle  and  SFP  data  cover  about  15%  of  the  35  year  orbit. 
Note  that  we  ignored  making  corrections  for  the  center  of  light 
motion  of  the  Ab  binary  because  the  largest  astrometric  shifts 
are  expected  to  be  small,  %0.5  mas. 

3.  VISUAL  ORBIT  OF  THE  CLOSE  SYSTEM 

The  interferometric  fringe  patterns  of  the  close  Abl,Ab2  pair 
of  stars  overlap  even  with  the  longest  baselines  available  at  the 
CHARA  Array,  so  we  cannot  use  the  SFP  method  to  measure 
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Table  2 

Visual  Orbital  Elements 


Element 

Aa,Ab  Orbit 

Abl,Ab2  Orbit 

P  (years)  . . . 

35.20  ±  1.45 

0.855333 

P  (days) . . . 

12855  ±  528 

312.40a 

T  (BY)  . . . 

1994.84  ±  1.69 

1996. 109a 

T  (HJD— 2.400,000)  . . . 

49662  ±616 

50123. 5a 

a  (mas) . . . 

54.5  ±  3.7 

3  Q+11 

i  (deg)  . . . 

46.2  ±  6.9 

51+17 

J1-51 

£2  (deg)  . . . 

255.2  ±  15.0 

25+3  b 
ZJ-35 

e  . . . 

0.489  ±0.081 

oa 

O)  (deg)  . . . 

70.4  ±  7.5 

180c 

Notes. 

a  Fixed  with  values  from  the  radial  velocity  orbit  (Table  6). 
b  Or  205*35  de§- 

c  Fixed  for  the  relative  orbit  of  Ab2  with  respect  to  Abl . 

the  relative  separation.  However,  the  interference  of  the  two 
fringe  patterns  of  the  inner  pair  causes  a  modulation  of  their 
combined  visibility  amplitude  with  changes  in  the  projected 
baseline  separation,  and  we  can  use  this  modulation  to  estimate 
the  binary  separation  projected  along  the  baseline  position  angle 
in  the  sky  (Hummel  et  al.  1998;  Boden  et  al.  2000;  Raghavan 
et  al.  2009).  The  calibration  of  visibility  is  aided  when  the  signal 
of  a  nearby  star  produces  an  SFP  that  can  be  used  to  calibrate 
the  visibility  of  the  central  binary.  The  details  of  this  method  are 
outlined  by  O’Brien  et  al.  (201 1). 

We  obtained  195  observations  of  HD  193322  over  24  nights 
between  2005  and  2010  using  the  CHARA  Array  Classic  beam 
combiner  (ten  Brummelaar  et  al.  2005).  The  observations  were 
made  with  the  near-IR  K'  filter  and  a  variety  of  different 
baselines.  These  measurements  were  a  series  of  approximately 
200  recorded  fringe  scans  sampled  at  a  frequency  of  150  Hz.  The 
scans  were  reduced  by  standard  techniques  (ten  Brummelaar 
et  al.  2005),  and  a  subset  of  100  scans  with  the  best  S/N 
were  selected.  We  fit  these  scans  with  fringe  patterns  for  both 
the  calibrator  (Aa)  and  target  (Ab)  (identified  according  to  the 
orbital  solution  from  Section  2  and  the  position  angle  of  the 
specific  baseline;  see  details  of  the  procedure  in  O’Brien  et  al. 
2011).  We  retained  only  those  visibility  measurements  for  which 
the  fractional  difference  between  a  first  estimate  and  final  mean 
was  less  than  20%,  and  these  best-case  values  were  used  to  form 
the  mean  visibility  for  each  component.  Note  that  we  used  the 
same  data  set  as  that  for  the  wide  pair  (Section  2),  but  due  to 
the  stronger  selection  criteria  not  all  data  sets  yielded  useful 
visibility  amplitudes.  We  next  determined  the  ratio  of  the  mean 
visibilities  of  the  target  and  calibrator.  This  observed  ratio  is 
related  to  the  ratio  of  the  individual  visibilities  for  the  target  and 
calibrator  by 

W  ^Aa  =  (FA,/FKh)  VAb,o/  VAa,o,  (1) 

where  f  Aa/  f  Ah  is  the  monochromatic  flux  ratio  in  the  K'  band. 
The  angular  diameter  of  the  calibrator  Aa  is  small  enough  that 
VAa  ^  1.0  for  our  observations,  but  we  estimated  the  single¬ 
star  visibilities  of  each  component  based  upon  the  projected 
baseline  of  observation  and  the  predicted  angular  diameters. 
A  value  for  the  flux  ratio  of  FAa/,FAb  =  ioJ°'4Amwide  =  0.92 
was  adopted  based  upon  a  fit  of  the  observed  ratios  and 
a  binary  model  (see  below),  and  this  parameter  essentially 
normalizes  the  target  visibility  so  that  the  upper  distribution 
of  the  visibilities  has  a  mean  of  one.  Note  that  we  expected 
that  the  ratio  would  be  f  Aa/  /*  Ab  Js  1  based  on  the  speckle  orbit 
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Figure  2.  Astrometric  orbit  of  the  Aa,Ab  pair  (in  units  of  arcsec).  The  dashed 
curve  represents  the  first  solution  (Hartkopf  et  al.  1993)  while  the  solid  curve  is 
the  new  solution  (Table  2).  The  dot-dashed  line  shows  the  line  of  the  nodes.  The 
filled  circles  represent  the  CHARA  Array  separated  fringe  packet  results  and 
open  circles  represent  the  speckle  observations.  Each  measurement  is  connected 
by  a  line  segment  to  the  calculated  position  for  the  time  of  observation.  Note 
that  north  is  down  and  east  is  to  the  right  in  this  figure,  and  the  directional  arc 
in  the  lower  right  comer  shows  the  counterclockwise  sense  of  orbital  motion. 

assignments  (Section  2),  but  we  suspect  that  this  difference  is 
probably  insignificant  given  the  uncertainties  in  the  component 
flux  fractions.  The  results  for  the  Ab  pair  are  given  in  Table  3 
which  lists  the  heliocentric  Julian  date  of  observation,  the 
corresponding  orbital  phase  in  the  Abl,Ab2  orbit  (Section  4), 
the  projected  baseline  Bp  and  position  angle  i jrp  of  observation, 
the  calibrated  visibility  and  its  associated  error,  and  the  observed 
minus  calculated  difference  O  —  C  in  visibility  from  the  adopted 
model  fit.  These  data  are  available  as  part  of  the  online  materials 
in  an  OIFITS  file. 

The  modulation  of  the  visibility  ratio  depends  on  the  known 
projected  baseline  length  and  position  angle  and  the  effec¬ 
tive  wavelength  of  the  K'  system,  plus  the  unknown  pro¬ 
jected  binary  separation  and  magnitude  differences,  A mciose  = 
—2.5  log(F’Ab2/ Fam)  and  A mwide  =  -2.5  log(FAa/FAb).  The 
latter  magnitude  difference  normalizes  the  visibility  according 
to  the  relation  given  above,  while  the  former  magnitude  differ¬ 
ence  sets  the  amplitude  of  the  visibility  modulation  with  the 
baseline  (Raghavan  et  al.  2009;  O’Brien  et  al.  201 1).  Following 
the  example  of  O’Brien  et  al.  (2011),  we  explored  the  orbital 
parameters  of  the  Abl,Ab2  pair  by  creating  a  set  of  model  visi¬ 
bilities  for  each  of  the  observed  times  and  baseline  parameters 
and  then  forming  the  x2  statistic  for  the  differences  between  the 
observed  and  model  visibilities.  The  solution  is  found  by  deter¬ 
mining  the  orbital  parameters  and  magnitude  differences  in  a 
high-resolution  grid  of  values  that  minimize  x2  -  For  this  appli¬ 
cation  of  the  method,  we  set  the  orbital  period  and  epoch  from 
the  spectroscopic  elements  for  the  circular  orbit  of  Abl,Ab2 
(Section  4)  and  then  made  a  grid  search  for  the  best-fit  values 
of  the  angular  semimajor  axis  a,  inclination  i,  and  longitude  of 
the  ascending  node  £2,  plus  the  magnitude  differences  Amc\ose 
and  A  m  wide- 
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Table  3  Table  3 

Visibility  Measurements  of  Abl,Ab2  (Continued) 


Date 

(HJD— 2,400,000) 

0 

(Close) 

B„ 

(m) 

fP 

(deg) 

V 

AV 

o-c 

Date 

(HJD— 2,400,000) 

<P 

(Close) 

BP 

(m) 

f p 
(deg) 

V 

AV 

O-C 

53591.721 

0.102 

98.8 

118.8 

1.035 

0.175 

0.045 

53914.829 

0.136 

97.6 

120.7 

1.029 

0.137 

0.036 

53591.756 

0.102 

104.3 

109.0 

0.968 

0.147 

-0.020 

53914.837 

0.136 

99.2 

118.1 

1.044 

0.144 

0.059 

53591.785 

0.102 

107.0 

101.9 

0.912 

0.125 

-0.024 

53914.842 

0.136 

99.9 

117.0 

0.940 

0.125 

-0.038 

53591.818 

0.102 

107.9 

93.9 

0.812 

0.108 

-0.039 

53914.850 

0.136 

101.4 

114.5 

0.925 

0.125 

-0.039 

53591.842 

0.102 

106.9 

88.5 

0.697 

0.095 

-0.104 

53914.859 

0.136 

102.7 

112.2 

0.937 

0.123 

-0.009 

53591.855 

0.102 

105.7 

85.3 

0.740 

0.097 

-0.038 

53914.865 

0.136 

103.6 

110.4 

0.872 

0.124 

-0.059 

53591.888 

0.102 

101.1 

77.2 

0.679 

0.093 

-0.073 

53950.685 

0.251 

88.2 

136.2 

0.896 

0.123 

0.077 

53591.912 

0.102 

96.2 

70.6 

0.667 

0.088 

-0.084 

53950.695 

0.251 

90.2 

132.7 

0.704 

0.093 

-0.095 

53638.741 

0.252 

170.6 

324.4 

0.893 

0.242 

0.137 

53950.701 

0.251 

91.5 

130.5 

0.773 

0.103 

-0.015 

53638.747 

0.252 

169.8 

323.5 

0.702 

0.119 

-0.057 

53950.708 

0.251 

92.9 

128.3 

0.770 

0.102 

-0.007 

53638.753 

0.252 

168.9 

322.4 

0.853 

0.153 

0.090 

53950.716 

0.251 

94.5 

125.7 

0.692 

0.090 

-0.075 

53638.757 

0.252 

168.2 

321.7 

0.797 

0.150 

0.031 

53950.724 

0.251 

96.1 

123.1 

0.759 

0.100 

-0.000 

53638.765 

0.252 

166.8 

320.4 

0.724 

0.129 

-0.047 

53950.730 

0.251 

97.3 

121.1 

0.687 

0.091 

-0.067 

53638.769 

0.252 

166.1 

319.8 

0.834 

0.160 

0.061 

53950.736 

0.251 

98.5 

119.3 

0.667 

0.089 

-0.085 

53638.775 

0.252 

164.9 

318.9 

0.802 

0.170 

0.025 

53950.744 

0.251 

100.0 

116.9 

0.733 

0.096 

-0.017 

53638.779 

0.252 

164.0 

318.3 

0.778 

0.173 

-0.001 

53950.755 

0.251 

101.7 

113.9 

0.760 

0.101 

0.010 

53638.782 

0.252 

163.3 

317.8 

0.940 

0.169 

0.159 

53950.764 

0.251 

103.0 

111.6 

0.678 

0.090 

-0.074 

53638.787 

0.252 

162.2 

317.1 

1.023 

0.145 

0.240 

53950.781 

0.251 

105.2 

107.1 

0.797 

0.105 

0.039 

53638.791 

0.252 

161.2 

316.5 

0.920 

0.174 

0.135 

53950.788 

0.251 

105.9 

105.3 

0.778 

0.104 

0.017 

53638.794 

0.252 

160.6 

316.1 

0.780 

0.162 

-0.007 

53950.803 

0.251 

107.1 

101.7 

0.849 

0.122 

0.082 

53638.798 

0.252 

159.4 

315.4 

0.884 

0.166 

0.096 

53950.814 

0.251 

107.7 

98.9 

0.826 

0.121 

0.055 

53638.804 

0.252 

158.0 

314.7 

0.817 

0.174 

0.027 

53950.819 

0.251 

107.8 

97.9 

0.729 

0.105 

-0.043 

53638.809 

0.252 

156.5 

314.0 

0.938 

0.155 

0.148 

53950.830 

0.251 

107.9 

95.3 

0.719 

0.103 

-0.054 

53638.813 

0.253 

155.4 

313.5 

0.931 

0.218 

0.140 

53950.838 

0.251 

107.8 

93.2 

0.774 

0.103 

-0.000 

53639.696 

0.255 

107.6 

91.8 

0.831 

0.118 

0.057 

53950.843 

0.251 

107.7 

92.2 

0.829 

0.114 

0.055 

53639.701 

0.255 

107.4 

90.6 

0.778 

0.119 

0.005 

53950.852 

0.251 

107.3 

90.1 

0.808 

0.119 

0.035 

53639.705 

0.255 

107.2 

89.6 

0.817 

0.120 

0.045 

53950.861 

0.251 

106.7 

88.0 

0.720 

0.096 

-0.051 

53639.713 

0.255 

106.7 

87.8 

0.901 

0.140 

0.130 

53950.866 

0.251 

106.2 

86.6 

0.839 

0.114 

0.069 

53639.717 

0.255 

106.3 

86.8 

0.775 

0.111 

0.005 

53950.875 

0.251 

105.3 

84.5 

0.785 

0.110 

0.018 

53639.723 

0.255 

105.7 

85.3 

0.765 

0.118 

-0.002 

53950.888 

0.251 

103.8 

81.5 

0.751 

0.104 

-0.012 

53639.729 

0.255 

105.0 

83.8 

0.594 

0.087 

-0.171 

53950.900 

0.251 

101.9 

78.4 

0.836 

0.112 

0.079 

53639.732 

0.255 

104.6 

83.1 

0.722 

0.122 

-0.041 

53950.904 

0.251 

101.3 

77.5 

0.708 

0.096 

-0.048 

53639.736 

0.255 

104.2 

82.2 

0.732 

0.106 

-0.031 

53982.762 

0.353 

174.5 

331.3 

1.133 

0.154 

0.199 

53639.739 

0.255 

103.8 

81.5 

0.829 

0.159 

0.068 

53982.769 

0.353 

174.0 

330.0 

1.068 

0.147 

0.135 

53639.743 

0.255 

103.2 

80.5 

0.785 

0.116 

0.026 

53982.777 

0.354 

173.2 

328.4 

1.070 

0.142 

0.137 

53639.746 

0.255 

102.6 

79.6 

0.888 

0.173 

0.130 

53982.782 

0.354 

172.7 

327.4 

1.088 

0.141 

0.156 

53639.750 

0.256 

102.0 

78.6 

0.834 

0.145 

0.078 

53982.791 

0.354 

171.8 

326.0 

1.017 

0.136 

0.083 

53639.754 

0.256 

101.4 

77.6 

0.638 

0.097 

-0.117 

53982.797 

0.354 

171.0 

324.9 

1.017 

0.135 

0.081 

53639.757 

0.256 

100.9 

76.9 

0.734 

0.110 

-0.020 

53982.805 

0.354 

169.8 

323.5 

1.008 

0.135 

0.069 

53639.767 

0.256 

99.0 

74.3 

0.792 

0.148 

0.040 

53982.814 

0.354 

168.5 

322.1 

0.975 

0.135 

0.031 

53639.769 

0.256 

98.5 

73.6 

0.729 

0.128 

-0.022 

53982.818 

0.354 

167.8 

321.3 

1.018 

0.136 

0.072 

53639.774 

0.256 

97.5 

72.3 

0.848 

0.198 

0.098 

53982.827 

0.354 

166.1 

319.9 

1.053 

0.138 

0.100 

53639.778 

0.256 

96.7 

71.3 

0.611 

0.116 

-0.139 

53982.849 

0.354 

161.4 

316.6 

0.677 

0.110 

-0.292 

53639.781 

0.256 

96.0 

70.3 

0.666 

0.113 

-0.085 

53982.854 

0.354 

160.2 

315.9 

1.053 

0.138 

0.080 

53639.785 

0.256 

95.1 

69.2 

0.832 

0.170 

0.081 

53982.863 

0.354 

157.7 

314.6 

1.040 

0.139 

0.060 

53639.788 

0.256 

94.4 

68.3 

0.917 

0.147 

0.166 

54273.899 

0.285 

105.6 

106.2 

0.747 

0.112 

-0.032 

53639.793 

0.256 

93.3 

66.9 

0.938 

0.198 

0.185 

54273.910 

0.285 

106.6 

103.5 

0.702 

0.103 

-0.079 

53639.797 

0.256 

92.3 

65.7 

1.270 

0.244 

0.515 

54273.914 

0.285 

106.9 

102.3 

0.795 

0.110 

0.014 

53639.800 

0.256 

91.7 

64.9 

1.379 

0.290 

0.622 

54273.930 

0.286 

107.7 

98.5 

0.798 

0.106 

0.018 

53639.808 

0.256 

89.7 

62.4 

1.003 

0.240 

0.240 

54273.938 

0.286 

107.9 

96.6 

0.789 

0.104 

0.011 

53639.812 

0.256 

88.6 

61.1 

0.776 

0.159 

0.008 

54273.948 

0.286 

107.9 

94.4 

0.789 

0.105 

0.013 

53639.816 

0.256 

87.4 

59.5 

0.922 

0.221 

0.148 

54273.958 

0.286 

107.7 

91.9 

0.879 

0.114 

0.107 

53639.822 

0.256 

86.0 

57.7 

1.040 

0.222 

0.258 

54273.963 

0.286 

107.5 

90.8 

0.794 

0.103 

0.024 

53893.877 

0.069 

315.5 

39.3 

1.154 

0.163 

0.357 

54273.979 

0.286 

106.3 

86.9 

0.847 

0.109 

0.085 

53893.886 

0.069 

318.1 

37.7 

0.989 

0.144 

0.209 

54273.987 

0.286 

105.6 

85.1 

0.785 

0.101 

0.027 

53893.889 

0.069 

319.0 

37.1 

0.893 

0.127 

0.117 

54273.997 

0.286 

104.4 

82.7 

0.719 

0.093 

-0.035 

53893.898 

0.069 

321.1 

35.5 

0.839 

0.130 

0.068 

54285.938 

0.324 

247.9 

8.0 

0.761 

0.169 

-0.130 

53893.912 

0.069 

323.8 

33.0 

1.062 

0.146 

0.287 

54288.939 

0.334 

247.9 

6.0 

1.072 

0.164 

0.090 

53893.915 

0.069 

324.3 

32.4 

1.088 

0.154 

0.309 

54288.986 

0.334 

248.0 

355.4 

0.860 

0.141 

-0.080 

53893.923 

0.069 

325.6 

30.8 

1.180 

0.160 

0.392 

54289.971 

0.337 

248.0 

358.2 

0.755 

0.115 

-0.195 

53893.926 

0.069 

325.9 

30.3 

1.031 

0.136 

0.240 

54289.975 

0.337 

248.0 

357.3 

0.700 

0.272 

-0.240 

53893.935 

0.069 

327.0 

28.6 

0.931 

0.130 

0.128 

54318.890 

0.429 

330.7 

2.5 

1.094 

0.144 

0.316 

53914.806 

0.136 

93.0 

128.0 

1.021 

0.134 

0.024 

54412.690 

0.730 

89.5 

62.2 

1.096 

0.226 

0.345 

53914.815 

0.136 

94.8 

125.1 

1.044 

0.160 

0.045 

54412.715 

0.730 

82.9 

53.5 

0.839 

0.141 

0.073 

53914.824 

0.136 

96.7 

122.1 

1.072 

0.141 

0.075 

54412.727 

0.730 

79.9 

49.1 

0.716 

0.115 

-0.065 
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Table  3 

(Continued) 


Date 

(HJD— 2,400,000) 

0 

(Close) 

BP 

(m) 

fP 

(deg) 

V 

AV 

O-C 

54412.739 

0.730 

76.8 

44.2 

0.801 

0.131 

-0.000 

54606.005 

0.348 

278.4 

143.2 

0.732 

0.191 

-0.218 

54657.944 

0.515 

267.3 

127.0 

0.731 

0.156 

-0.092 

54657.959 

0.515 

262.1 

124.1 

0.986 

0.191 

0.074 

54657.968 

0.515 

258.6 

122.5 

1.115 

0.181 

0.161 

54692.830 

0.626 

272.2 

130.7 

0.684 

0.096 

-0.147 

54692.837 

0.626 

270.6 

129.3 

0.727 

0.096 

-0.143 

54692.889 

0.627 

330.7 

177.4 

1.008 

0.146 

0.033 

54692.897 

0.627 

330.7 

175.6 

1.068 

0.144 

0.099 

54692.905 

0.627 

330.7 

173.9 

0.761 

0.112 

-0.196 

54692.912 

0.627 

330.7 

172.0 

0.702 

0.127 

-0.232 

54692.946 

0.627 

330.4 

164.3 

0.702 

0.196 

-0.087 

54692.960 

0.627 

330.1 

161.2 

0.580 

0.128 

-0.174 

54759.629 

0.840 

275.2 

134.0 

0.983 

0.168 

0.063 

54759.667 

0.840 

330.7 

186.2 

0.997 

0.175 

0.193 

54759.677 

0.840 

330.7 

183.9 

1.271 

0.183 

0.503 

54759.687 

0.840 

330.7 

181.6 

1.092 

0.203 

0.340 

54759.696 

0.840 

330.7 

179.5 

1.380 

0.215 

0.623 

54759.728 

0.841 

238.0 

115.8 

1.088 

0.206 

0.332 

54759.765 

0.841 

330.3 

163.6 

0.809 

0.133 

-0.151 

54759.790 

0.841 

329.5 

158.0 

0.816 

0.133 

-0.160 

54983.996 

0.558 

277.2 

137.8 

0.640 

0.091 

-0.172 

54984.002 

0.558 

276.7 

136.6 

0.766 

0.100 

-0.076 

54984.865 

0.561 

262.4 

100.2 

1.204 

0.160 

0.239 

54984.871 

0.561 

267.1 

98.6 

1.186 

0.156 

0.260 

54984.876 

0.561 

270.9 

97.3 

1.175 

0.163 

0.288 

54984.882 

0.561 

274.5 

96.1 

0.926 

0.124 

0.079 

54984.887 

0.561 

278.1 

94.8 

0.926 

0.121 

0.122 

54984.892 

0.561 

281.4 

93.6 

0.785 

0.105 

0.013 

54984.904 

0.561 

276.2 

159.5 

1.145 

0.151 

0.179 

54984.907 

0.561 

276.5 

158.2 

1.150 

0.151 

0.194 

54984.916 

0.561 

276.9 

156.2 

1.012 

0.139 

0.077 

54984.921 
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Figure  3.  Gray-scale  representation  of  the  minimum  /  2  over  the  full  range  in  Q 
for  visibility  models  of  the  Abl  ,Ab2  binary  as  a  function  of  orbital  inclination  i 
and  angular  semimajor  axis  a.  Within  this  numerical  grid  of  (/,  a ),  the  minimum 
is  x2  =  167  (black)  and  the  maximum  is  /2  =  529  (white)  for  a  sample  of 
195  measurements  and  five  fitting  parameters.  The  solid  lines  indicate  loci  of 
constant  M(Abl)  (10,  20,  30,  40,  and  50  Mq  from  bottom  to  top),  while  the 
dashed  lines  represent  loci  of  constant  M( Ab2)  (5,  10,  15,  and  20  Mq  from 
bottom  to  top),  all  for  an  assumed  distance  of  741  pc.  The  location  of  the 
adopted  solution  is  marked  by  a  plus  sign. 


We  found  that  the  solutions  always  arrived  at  similar  estimates 
for  the  magnitude  differences,  so  we  set  these  magnitude 
differences  and  performed  a  grid  search  over  i  and  a ,  the  two 
parameters  of  physical  interest.  For  each  selection  of  (;,  a), 
we  determined  the  best  fit  for  the  sky  orientation  parameter  £2 
over  the  full  range  of  values  in  steps  of  2°.  The  resulting  y2 
estimates  are  plotted  in  a  gray-scale  diagram  in  Figure  3  in  the 
(;,  a)  plane  for  grid  increments  of  A;  =  2°  and  A  a  =  0.05  mas. 
Here  intensity  is  scaled  between  the  lowest  (black)  and  highest 
(white)  y2  over  the  grid.  If  we  assume  a  distance  d  from  the 
cluster  fitting  results  of  Roberts  et  al.  (2010),  then  Kepler’s  third 
law  relates  the  known  period  P,  the  total  mass,  and  a  by 


(M(Ab  1 )  +  M(  Ab2))/M0 


(ad)2  /  a  \3/  d  \3 

P2  vl.22mas/  \741pc/ 


(2) 

Next,  we  can  use  the  spectroscopic  semiamplitude  K  for  com¬ 
ponent  Abl  (Section  4,  Table  6)  to  derive  a  relation  for  the  mass 
of  Ab2  as  a  function  of  i  and  a. 


M(  Ab2)  = 


a2  K(  1  -  e2)l/2 
Psini  29.8  km  s-1 


0.455  Mq - 

sin  i 


Then  we  can  find  the  mass  of  Abl  from  a  relation  for  the  mass 
ratio, 


M(Abl)  29.8  km  s  1  adsini  ( 

M (Ab2)  "  K  ~ 

asm;  /  d  \ 

0.819 mas  \741  pc/  ^  ^ 

Thus,  each  point  in  the  (;,  a)  plane  is  associated  with  specific 
masses  M( Abl)  and  M(Ab2),  and  we  can  use  the  relations 
above  to  construct  loci  of  constant  primary  and  secondary  mass 
in  Figure  3  (shown  by  solid  and  dashed  lines,  respectively). 

Inspection  of  Figure  3  indicates  that  there  are  two  broad 
valleys  in  the  (;,  a)  plane  where  the  fits  are  relatively  good,  one 
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Figure  4.  Plots  of  the  calibrated  (plus  sign)  and  model  (diamond)  interferometric  visibilities  for  component  Abl  ,Ab2  for  each  of  the  first  12  nights  of  observation  with 
the  CHARA  Array. 


with  i  <  90°  for  counterclockwise  motion  in  the  sky  and  another 
with  i  >  90°  for  clockwise  motion.  Within  these  valleys  there 
are  three  locations  with  comparable  minima,  but  all  of  these 
are  associated  with  extreme  masses:  (/,  a)  —  (66°,  4.6  mas) 
(X2  =  210)  and  (/,  a)  =  (1 18°,  4.7  mas)  (x2  =  215)  where 
the  masses  are  too  high  and  ( i ,  a)  —  (38°,  3.3  mas)  (x2  =  167) 
where  the  masses  are  too  low  (see  Section  6  below).  We  think 
that  the  i  <  90°  valley  probably  represents  the  best  family  of 
solutions,  since  the  trends  in  x2  are  more  or  less  continuous 
there  as  expected.  Between  a  —  2.7  and  5.0  mas,  the  valley 
floor  never  rises  above  a  reduced  chi-square  of  x2  =  1.18  (with 
190  degrees  of  freedom,  equal  to  195  measurements  minus  five 
fitting  parameters).  Although  a  purely  statistical  assessment 
would  restrict  the  solution  space  to  the  valley  region  around 
(/,  a)  —  (38°,  3.3  mas),  the  fact  that  the  reduced  chi-square  is 
close  to  unity  along  the  length  of  the  valley  suggests  that  at 
this  stage  it  is  premature  to  rule  out  any  of  this  solution  space. 
In  Section  6  below,  we  present  several  lines  of  argument  that 
indicate  that  the  actual  solution  lies  in  the  mid-range  of  this 
valley  at  a  —  3.85  mas  (x2  =  225),  so  we  will  tentatively  adopt 
this  value  and  present  the  associated  solution  for  the  other  orbital 
parameters  in  Table  2,  Column  3.  The  errors  associated  with  /, 
a,  and  £2  reported  in  Table  2  correspond  to  their  range  over  the 
length  of  the  valley  from  a  =  2. 7-5.0  mas.  Note  that  because 
the  visibility  oscillation  depends  on  the  absolute  value  of  the 


projected  separation,  there  is  a  180°  ambiguity  in  our  derived 
value  of  £2.  We  found  that  the  best-fit  magnitude  differences  are 
Amciose  =  2. 1 1  ±  0.06  mag  and  Amwl(ic  =  0.086  ±  0.012  mag 
(Ab  brighter  than  Aa).  In  order  to  show  how  well  the  model  and 
observed  visibilities  agree,  we  plot  the  individual  and  calculated 
visibilities  for  this  solution  for  each  night  in  Figures  4  and  5, 
and  we  find  that  the  fits  are  satisfactory  for  most  of  the  nights. 

The  projection  of  the  orbit  on  the  sky  for  this  solution  is 
illustrated  in  Figure  6  where  filled  circles  indicate  the  calculated 
positions  at  the  times  of  observation.  The  distribution  of  the 
observations  in  orbital  phase  appears  to  constrain  the  minor 
axis  of  the  projected  ellipse  better  than  the  major  axis.  With  the 
minor  axis  fixed,  the  major  axis  will  vary  with  inclination  as 
a(minor)  ~  a  cos  /  or  a  oc  sec  /,  and  this  relation  approximately 
describes  the  position  of  the  x2  valley  in  Figure  3.  The  orbital 
orientations  of  the  wide  and  close  orbits  appear  quite  different 
(compare  Figures  2  and  6),  and  the  mutual  inclination  of  the 
orbital  planes  0  is  given  by  Fekel  (1981)  as 

COS  (j)  —  COS  /close  COS  /wide  T  sill  /close  sin  /wide 

X  COS(£2close  £2wide )■  (5) 

where  the  second  term  may  assume  either  sign  because  of 
the  180°  ambiguity  in  the  determination  of  £2ci0se-  The  two 
solutions,  <p  =  38°  and  85°,  indicate  that  the  two  orbits  are 
probably  far  from  coplanar  (<f>  =  0°). 
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Figure  5.  Plots  of  the  calibrated  (plus  sign)  and  model  (diamond)  interferometric  visibilities  for  component  Abl,Ab2  for  each  of  the  last  12  nights  of  observation  with 
the  CHARA  Array.  Note  that  the  self-calibration  method  used  for  these  data  is  extremely  seeing  dependent  and  this  can  cause  large  differences  between  the  model 
and  the  data  on  some  evenings.  This  is  not  unusual  in  interferometric  data  of  low  S/N. 


ARA  (mas) 

Figure  6.  Astrometric  orbit  of  the  Ab  1 .  Ab2  pair  in  the  same  format  as  Figure  2 
(but  now  in  units  of  milliarcsecond  (mas))  based  upon  the  CHARA  Array 
visibility  measurements.  Filled  circles  indicate  the  calculated  positions  at  the 
times  of  observation. 


4.  RADIAL  VELOCITIES  OF  AB1 

Spectroscopy  potentially  offers  us  the  means  to  determine 
the  masses  and  spectral  properties  of  the  components.  However, 
because  the  stars  are  so  close,  conventional,  ground-based 
spectroscopy  records  the  flux  of  all  three  stars  (usually  plus 
component  B  at  a  separation  of  2" 7 ;  Turner  et  al.  2008),  and  since 
their  orbital  Doppler  shifts  are  comparable  to  the  line  widths,  the 
resulting  line  blending  problem  is  daunting.  Nevertheless,  the 
spectral  properties  of  the  components  are  sufficiently  different  in 
this  case  that  we  may  attempt  radial  velocity  measurements.  The 
appearance  of  the  optical  spectrum  is  dominated  by  a  narrow- 
lined  component  that  corresponds  to  the  primary  star  in  the 
close  orbit,  Abl  (McKibben  et  al.  1998).  The  secondary  in  the 
close  orbit,  Ab2,  is  fainter  and  contributes  little  to  the  composite 
spectrum  (Section  5).  Furthermore,  there  is  a  very  broad  lined 
component  that  appears  to  follow  the  motion  of  Aa  in  the  wide 
orbit  (Section  5).  Because  the  lines  of  Aa  are  so  broad  and 
shallow,  they  essentially  act  to  depress  the  continuum  in  the 
vicinity  of  the  narrow  lines  of  Abl,  and  since  the  velocity  range 
of  Ab  1  is  smaller  than  the  full  width  of  the  lines  of  Aa  in  general, 
the  presence  of  the  broad  component  has  little  influence  on 
velocity  measurements  of  Abl  (but  see  a  discussion  of  blending 
effects  below).  Here  we  present  radial  velocities  for  the  Abl 
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Source 

Number 

Date 

(BY) 

Range 

(A) 

Resolving  Power 
(X/AX) 

N 

Observatory /Telescope  / 
Spectrograph 

i... 

1995.5 

58445904 

26000 

1 

Ritter/ 1  m/Echellea 

2... 

1995.6 

5572-5895 

22200 

1 

KPNO/0.9  m/Coude 

3... 

1995.6 

6434-6751 

31000 

1 

KPNO/O.9  m/Coude 

4... 

1998.7 

6314-6978 

12200 

1 

KPNO/0.9  m/Coude 

5... 

1999.8 

5401-6735 

5600 

4 

KPNO/O.9  m/Coude 

6... 

2000.7 

6443-7108 

12500 

1 

KPNO/0.9  m/Coude 

7... 

2001.0 

6443-7108 

12500 

3 

KPNO/O.9  m/Coudeb 

8... 

2002.4 

4692-6018 

4900 

1 

KPNO/O.9  m/Coude 

9... 

2002.4 

5980-7313 

6100 

1 

KPNO/0.9  m/Coude 

10... 

2004.7 

4000-6800 

34200 

2 

OHP/1.9  m/Elodiec 

11... 

2004.8 

6466-7176 

7900 

1 

KPNO/O.9  m/Coude 

12... 

2005.9 

4236-4587 

10300 

2 

KPNO/O.9  m/Coude 

13... 

2006.8 

6466-7176 

7900 

2 

KPNO/0.9  m/Coude 

14... 

2006.8 

4236-4587 

10300 

2 

KPNO/0.9  m/Coude 

15... 

2008.6 

4465~4586 

76300 

2 

CFHT/3.6  m/ESPaDOnS 

16... 

2009.9 

4000-4720 

75900 

1 

NOT/2.6  m/FIES 

17... 

2010.5 

3994-4663 

5700 

2 

KPNO/4.0  m/R-C 

18... 

2010.5 

3873^4540 

6400 

1 

Lowell/ 1.8  m/DeVeny 

19... 

2010.6 

4292^4670 

4300 

1 

DAO/1.8  m/Cassegrain 

20... 

2010.6 

3873^4540 

6400 

1 

Lowell/1.8  m/DeVeny 

Notes. 

a  http://astrol.panet.utoledo.edu/~wwritter/archive/PREST-archive.html 
b  http://www.noao.edu/cllib/ 
c  http://atlas.obs-hp.fr/elodie/ 


component  and  show  that  they  represent  the  sum  of  orbital 
motions  in  both  the  close  and  wide  systems. 

We  collected  3 1  new  spectra  for  measurement  from  sources 
which  are  summarized  in  Table  4.  The  columns  list  a  source 
number  (for  identification  with  the  specific  radial  velocities 
listed  in  Table  5),  date  of  observation(s),  spectral  range  used  in 
the  measurement,  the  spectral  resolving  power,  number  of  spec¬ 
tra  made  at  that  time,  and  the  observatory,  telescope  aperture, 
and  spectrograph  of  origin.  We  obtained  most  of  these  spectra 
in  runs  at  the  Kitt  Peak  National  Observatory  0.9  m  coude  feed 
and  4  m  Mayall  Telescopes,  the  3.6  m  Canada-France-Hawaii 
Telescope  (CFHT),  the  2.5  m  Nordic  Optical  Telescope,  the 
Lowell  Observatory  1.8  m  Perkins  Telescope,  and  the  Herzberg 
Institute  of  Astrophysics,  Dominion  Astrophysical  Observatory 
1.8  m  Plaskett  Telescope.  These  were  augmented  with  publicly 
available  spectra  from  the  archives  of  the  University  of  Toledo 
Ritter  Observatory  1.0  m  telescope  (Morrison  et  al.  1997),  the 
Observatoire  de  Haute-Provence  1.9  m  telescope  and  ELODIE 
spectrograph  (Moultaka  et  al.  2004),  and  the  Indo-U.S.  Library 
of  Coude  Feed  Stellar  Spectra  (Valdes  et  al.  2004;  made  with 
the  KPNO  0.9  m  coude  feed  telescope).  All  these  spectra  were 
reduced  by  standard  techniques  and  transformed  to  a  continuum 
normalized  flux  representation  on  a  heliocentric,  log  wavelength 
grid.  Atmospheric  telluric  lines  were  removed  from  the  red  spec¬ 
tra  by  division  with  a  pure  atmospheric  spectrum.  This  was  done 
by  creating  a  library  of  spectra  from  each  run  of  a  rapidly  ro¬ 
tating  A  star  (usually  £  Aql),  removing  the  broad  stellar  fea¬ 
tures  from  these,  and  then  dividing  each  target  spectrum  by  the 
modified  atmospheric  spectrum  that  most  closely  matched  the 
target  spectrum  in  a  selected  region  dominated  by  atmospheric 
absorptions. 

The  spectra  form  a  diverse  collection  with  a  wide  range  in 
resolving  power  and  wavelength  coverage.  In  order  to  measure 
radial  velocities  in  a  consistent  way,  we  cross-correlated  each 


of  the  spectra  with  a  standard,  model  spectrum  (rest  frame) 
from  the  grid  of  synthetic  spectra  from  Lanz  &  Hubeny  (2003). 
From  an  initial  inspection  of  the  observations,  we  selected  a 
model  with  Galactic  abundances,  an  effective  temperature  of 
Teff  =  34.8  kK,  a  gravity  log  of  g  =  4.0,  a  projected  rotational 
velocity  of  V  sin  i  =  50  km  s_1,  a  wavelength-dependent  limb- 
darkening  coefficient  from  Wade  &  Rucinski  (1985),  and  an 
instrumental  broadening  appropriate  for  the  specific  observed 
spectrum.  The  cross-correlations  were  generally  made  over  the 
wavelength  range  given  in  Table  4,  although  in  some  cases 
regions  with  strong  interstellar  features  were  omitted.  The 
resulting  cross-correlation  functions  were  always  single-peaked, 
and  we  measured  the  radial  velocity  and  its  associated  error 
using  the  method  of  Zucker  (2003).  The  results  are  presented  in 
Table  5  which  lists  the  heliocentric  Julian  date  of  mid-exposure, 
the  corresponding  orbital  phases  in  the  close  and  wide  systems 
(see  below),  the  measured  radial  velocity  and  its  associated 
error,  a  correction  term  for  line  blending  effects,  the  observed 
minus  calculated  velocity  residual  from  the  fit  (see  below),  and 
the  observation  source  number  from  Table  4.  Note  that  for 
completeness  we  have  included  in  Table  5  velocities  published 
earlier  by  McKibben  et  al.  (1998;  indicated  by  a  0  in  the  final 
column). 

Although  the  line  blending  effects  from  the  spectral  compo¬ 
nents  of  the  other  stars  are  generally  small,  they  tend  to  bias 
the  measurements  toward  the  systemic  velocity  and  lead  to  a 
slight  underestimate  of  the  orbital  semiamplitude.  The  radial 
velocity  offset  caused  by  line  blending  will  depend  on  the  char¬ 
acter  and  velocity  shift  of  each  component,  the  spectral  features 
measured,  and  the  spectral  resolution  of  the  observation.  In  or¬ 
der  to  make  a  simple  correction  for  line  blending  effects  we 
adopted  the  following  procedure  for  each  observation.  We  first 
determined  model  synthetic  spectra  for  each  stellar  component 
(Section  5)  for  the  spectral  range  and  instrumental  broadening 
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Table  5 

Radial  Velocity  Measurements  for  Abl 


Date 

(HJD— 2,400,000) 

0 

(Close) 

<P 

(Wide) 

vr 

(km  s'1) 

a(Vr) 

(km  s'1) 

A  Vr  (Blend) 
(km  s'1) 

0  -  C 
(km  s'1) 

Source 

Number3 

22815.476 

0.588 

0.966 

-16.1 

8.0 

-2.3 

-4.6 

0 

23952.434 

0.227 

0.054 

3.1 

8.0 

1.7 

10.4 

0 

23962.408 

0.259 

0.055 

3.8 

8.0 

1.5 

15.2 

0 

24254.017 

0.192 

0.078 

8.8 

8.0 

2.3 

14.1 

0 

24424.648 

0.739 

0.091 

-17.1 

8.0 

-2.3 

-6.8 

0 

24668.949b 

0.521 

0.110 

1.1 

8.0 

-0.5 

33.2 

0 

24673. 884b 

0.536 

0.110 

-6.3 

8.0 

-2.7 

23.3 

0 

24675.927 

0.543 

0.111 

-17.0 

8.0 

-3.1 

11.9 

0 

39002.424 

0.402 

0.225 

-26.0 

8.0 

-15.9 

-13.5 

0 

40044.910 

0.739 

0.306 

-14.9 

1.4 

-1.0 

-4.6 

0 

40065.488 

0.805 

0.308 

0.0 

8.0 

3.1 

5.7 

0 

40347.870 

0.708 

0.330 

-11.0 

1.1 

-1.0 

2.9 

0 

43741.479 

0.571 

0.594 

-21.5 

1.3 

-0.0 

1.7 

0 

43741.516 

0.571 

0.594 

-21.7 

1.3 

-0.0 

1.5 

0 

43771.710 

0.668 

0.596 

-17.7 

2.1 

-1.5 

-4.7 

0 

43772.750 

0.671 

0.596 

-16.4 

3.1 

-1.4 

-3.7 

0 

43777.740 

0.687 

0.596 

-8.1 

3.1 

-0.4 

3.7 

0 

44051.110 

0.562 

0.618 

-21.9 

1.3 

-0.0 

1.2 

0 

44087.791 

0.680 

0.621 

-11.8 

1.3 

0.0 

0.8 

0 

44593.740 

0.299 

0.660 

-4.5 

1.3 

0.4 

4.9 

0 

45659.984 

0.712 

0.743 

-0.2 

1.3 

0.3 

5.3 

0 

45991.853 

0.775 

0.769 

8.6 

1.3 

0.6 

5.4 

0 

46606.058 

0.741 

0.816 

-3.2 

3.3 

0.5 

-3.5 

0 

46607.015 

0.744 

0.816 

-1.9 

2.2 

0.6 

-2.6 

0 

46608.051 

0.747 

0.817 

-2.4 

2.1 

0.6 

-3.6 

0 

46609.062 

0.750 

0.817 

-0.9 

1.1 

0.6 

-2.4 

0 

46612.036 

0.760 

0.817 

3.1 

1.4 

0.8 

0.5 

0 

46985.281 

0.955 

0.846 

25.0 

1.3 

0.4 

2.0 

0 

46986.266 

0.958 

0.846 

21.5 

1.3 

0.4 

-1.6 

0 

46986.660 

0.959 

0.846 

23.2 

1.3 

0.4 

0.0 

0 

46986.707 

0.959 

0.846 

21.1 

1.3 

0.4 

-2.1 

0 

46986.778 

0.959 

0.846 

22.4 

1.3 

0.4 

-0.8 

0 

46988.682 

0.966 

0.846 

23.3 

1.3 

0.4 

-0.1 

0 

46988.724 

0.966 

0.846 

24.2 

1.3 

0.4 

0.8 

0 

46988.769 

0.966 

0.846 

21.6 

1.3 

0.4 

-1.8 

0 

46988.815 

0.966 

0.846 

22.6 

1.3 

0.4 

-0.8 

0 

47773.924 

0.479 

0.907 

-14.1 

1.3 

-0.0 

1.6 

0 

47773.969 

0.479 

0.907 

-15.0 

1.3 

-0.0 

0.7 

0 

49231.714 

0.145 

0.021 

10.9 

2.9 

5.2 

7.1 

0 

49236.773 

0.162 

0.021 

0.0 

2.1 

2.0 

-5.1 

0 

49614.855 

0.372 

0.050 

-24.7 

5.4 

0.5 

-1.4 

0 

49840.255 

0.093 

0.068 

8.5 

1.3 

0.2 

0.9 

0 

49842.803 

0.102 

0.068 

5.8 

4.0 

8.0 

6.7 

0 

49843.780 

0.105 

0.068 

10.4 

4.2 

10.1 

13.6 

0 

49916.748 

0.338 

0.074 

-23.1 

5.1 

-0.4 

-2.1 

1 

49942.734 

0.421 

0.076 

-30.0 

0.8 

-0.6 

-1.6 

2 

49942.886 

0.422 

0.076 

-23.9 

0.9 

-6.8 

-1.6 

3 

49982.849 

0.550 

0.079 

-28.0 

1.3 

0.2 

2.8 

0 

49985.625 

0.559 

0.079 

-16.1 

2.9 

-18.7 

-4.5 

0 

50059.437b 

0.795 

0.085 

16.3 

6.5 

11.7 

33.0 

0 

51056.748 

0.987 

0.163 

10.3 

1.9 

5.0 

6.1 

4 

51466.709 

0.300 

0.195 

-17.3 

2.0 

-1.4 

-0.6 

5 

51467.788 

0.303 

0.195 

-12.3 

2.0 

-1.2 

5.0 

5 

51467.795 

0.303 

0.195 

-12.1 

2.0 

-1.2 

5.2 

5 

51468.754 

0.306 

0.195 

-17.8 

2.2 

-1.5 

-0.4 

5 

51817.657 

0.423 

0.222 

-29.3 

1.4 

-3.4 

-2.8 

6 

51888.616 

0.650 

0.227 

-16.8 

1.2 

-1.7 

5.1 

7 

51893.557 

0.666 

0.228 

-14.1 

2.5 

-1.3 

6.4 

7 

51895.593 

0.672 

0.228 

-11.7 

2.9 

-1.0 

8.3 

7 

52430.950b 

0.386 

0.270 

-11.3 

2.8 

-2.2 

12.8 

8 

52436.914 

0.405 

0.270 

-26.8 

9.5 

-5.7 

-4.5 

9 

53246.460 

0.997 

0.333 

8.2 

0.9 

1.3 

-2.3 

10 

53247.481 

1.000 

0.333 

8.4 

1.1 

1.4 

-2.0 

10 

53290.656 

0.138 

0.336 

-1.8 

2.3 

0.9 

-5.3 

11 

53683.603 

0.396 

0.367 

-21.4 

1.0 

-2.4 

1.6 

12 
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Date 

<t> 

<t> 

Vr  Cr(Vr) 

A  Fr  (Blend) 

o-c 

1  Source 

(HJD— 2,400,000) 

(Close) 

(Wide) 

(km  s'1)  (km  s'1) 

(km  s'1) 

(km  s~ 

[)  Number2 

53684.593 

0.399 

0.367 

-21.3 

1.0 

-2.4 

2.1 

12 

54019.652 

0.472 

0.393 

-30.1 

2.3 

-3.4 

-4.4 

13 

54024.715 

0.488 

0.394 

-26.9 

2.3 

-3.1 

-0.7 

13 

54029.707 

0.504 

0.394 

-25.9 

0.9 

-3.1 

0.3 

14 

54031.627 

0.510 

0.394 

-23.9 

1.0 

-2.9 

2.5 

14 

54675.892 

0.572 

0.444 

-27.3 

0.7 

-2.2 

-3.2 

15 

54675.917 

0.572 

0.444 

-27.3 

0.6 

-2.2 

-3.2 

15 

55146.400 

0.078 

0.481 

9.8 

0.2 

2.2 

-0.2 

16 

55366.903 

0.784 

0.498 

7.3 

2.5 

2.6 

11.7 

17 

55369.913 

0.794 

0.498 

2.6 

2.5 

1.9 

4.9 

17 

55383.939 

0.839 

0.499 

6.9 

1.3 

3.4 

5.3 

18 

55402.849 

0.899 

0.501 

4.2 

2.4 

5.3 

-1.3 

19 

55402.871 

0.899 

0.501 

3.1 

1.4 

2.9 

-4.9 

20 

Notes. 

3  0:  McKibben  et  al.  (1998);  1-20:  see  Table  4. 
b  Assigned  zero  weight  in  the  orbital  solution. 

Table  6 

Radial  Velocity  Orbital  Elements  for  Abl 

Orbital 

Aa,Ab  System 

Aa,Ab  System 

Abl,Ab2  System 

Abl,Ab2  System 

Element 

(No  Correction) 

(Blend  Correction) 

(No  Correction) 

(Blend  Correction) 

P  (years)  . . . 

35.20a 

35.20a 

0.85543  ±  0.00030 

0.85533  ±  0.00029 

P  (days) . . . 

12855a 

12855a 

312.44  ±0.11 

312.40  ±0.10 

T  (BY)  . . . 

1993.1  ±0.4 

1992.9  ±  0.3 

1996.111  ±  0.005 

1996.109  ±0.004 

T  (HJD— 2,400.000) . . 

49030  ±  148 

48966  ±  103 

50124.1  ±  1.7 

50123.5  ±  1.5 

K  (km  s'1)  ... 

8.6  ±0.6 

8.7  ±  0.4 

19.2  ±  0.4 

21.1  ±0.4 

y  (km  s-1) . . . 

-4.3  ±0.4 

-4.7  ±  0.4 

-4.3  ±  0.4 

-4.7  ±  0.4 

e  . . . 

0.489a 

0.489a 

0.0 

0.0 

co  (deg)  . . . 

70.4a 

70. 4a 

f(m)  (  Mq)  . . . 

0.56  ±0.12 

0.58  ±  0.08 

0.230  ±  0.015 

0.306  ±  0.019 

a\  sini  (106  km) . . . 

1327  ±  94 

1341  ±  65 

82.5  ±  1.8 

90.8  ±  1.9 

rms  (km  s-1) . . . 

3.0 

3.1 

3.0 

3.1 

Note.  a  Fixed  with  values  from  the  visual  orbit  (Table  2). 


of  the  observation.  The  models  for  components  Aa,  Ab2,  and 
B  were  co-added  according  to  the  adopted  fluxes  and  to  the 
Doppler  shifts  for  the  time  of  observation.  Then  we  formed  a 
series  of  model  spectra  by  adding  in  component  Ab  1  for  a  grid 
of  assigned  velocity  offsets,  and  we  measured  the  radial  veloc¬ 
ity  in  these  composite  spectra  using  the  same  cross-correlation 
method  applied  to  the  observations.  This  led  to  a  relation  be¬ 
tween  the  actual  and  measured  radial  velocity  for  each  obser¬ 
vation,  and  we  interpolated  within  this  relation  at  the  observed 
radial  velocity  to  determine  the  offset  correction  for  blending 
AV, (blend)  =  V).  (actual)  —  \/,( measured),  which  is  given  in 
Table  5,  Column  6.  The  average  of  the  absolute  value  of  the 
offset  correction  for  blending  is  small,  2.3  km  s'1,  but  the  in¬ 
dividual  offset  corrections  are  larger  for  the  lower  resolution 
spectra  where  line  blending  is  more  severe. 

The  velocities  of  component  Ab  1  depend  on  its  orbital  motion 
in  the  close  binary  plus  the  motion  of  the  Abl,Ab2  center  of 
mass  in  the  orbit  of  the  Aa,Ab  system.  Our  first  solutions  for 
the  orbital  motion  of  the  close  binary  clearly  showed  long-term 
variations  in  the  residuals  that  followed  the  motion  predicted  for 
Ab  in  the  wide  orbit.  Thus,  we  fit  the  observed  radial  velocity 
variations  as  the  sum  of  motions  in  the  close  and  wide  binaries. 
This  was  done  iteratively  using  the  orbital  fitting  program  of 
Morbey  &  Brosterhus  (1974).  We  first  made  a  general  fit  of  the 


velocities  for  the  close  system,  and  then  we  made  a  constrained 
fit  of  the  velocity  residuals,  by  fixing  P,  e,  and  co  from  the 
visual  orbit  of  Aa,Ab  to  find  the  semiamplitude  K  and  epoch 
of  periastron  T  for  the  wide  system.  The  resulting  solution  of 
the  long-period  orbit  was  then  used  to  correct  the  observed 
velocities  for  motion  in  the  wide  orbit,  and  a  new  solution  was 
found  for  the  close  orbit.  This  procedure  quickly  converged 
to  yield  the  orbital  elements  given  in  Table  6.  Note  that  we 
assigned  each  measurement  a  weight  proportional  to  <j(Vr)~ 2  in 
making  the  fits,  and  we  zero  weighted  four  measurements  that 
had  unusually  large  residuals  from  the  final  fit  (dates  indicated 
in  Table  5  and  shown  as  open  circles  in  Figure  7).  Table  6  lists 
the  solutions  both  with  and  without  application  of  the  offset 
correction  for  line  blending,  and  they  are  generally  very  similar 
except  for  the  slightly  larger  semiamplitude  K  that  results  when 
accounting  for  line  blending.  Since  the  line  blending  problem 
is  significant,  we  adopt  the  corrected  velocity  solutions  that  are 
illustrated  in  Figure  7  (close  orbit)  and  Figure  8  (wide  orbit)  and 
that  form  the  basis  for  the  residuals  O  —  C  given  in  Table  5, 
Column  7.  We  found  that  the  eccentricity  associated  with  the 
close  orbit  is  not  statistically  different  from  zero  according  to 
the  criterion  of  Lucy  &  Sweeney  (1971),  so  we  present  circular 
elements  in  Table  6  (where  the  epoch  T  is  defined  as  the  time  of 
maximum  radial  velocity  or,  equivalently,  the  time  of  crossing 
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Figure  7.  Derived  radial  velocity  curve  of  Abl  (solid  line)  in  the  312  day  orbit. 
Open  circles  indicate  those  four  measurements  assigned  zero  weight  in  the 
solution.  Phase  zero  corresponds  to  the  time  of  maximum  radial  velocity  (star 
crossing  the  ascending  node)  in  this  circular  orbit. 


ORBITAL  PHASE 

Figure  8.  Derived  radial  velocity  curve  of  Abl,Ab2  (low  amplitude,  solid  line) 
in  the  34  year  orbit.  Open  circles  indicate  measurements  of  the  broad-lined 
Aa  component  and  its  preliminary  radial  velocity  curve  (large  amplitude,  solid 
line).  Phase  zero  corresponds  to  the  time  of  periastron. 


the  ascending  node).  The  long  orbital  period  of  the  wide  system, 
P  =  35  years,  places  HD  193322  among  the  top  1%  of  known 
spectroscopic  binaries  with  very  long  periods  (Pourbaix  et  al. 
2004). 

5.  SPECTROSCOPIC  PROPERTIES 

The  two  brightest  components  of  HD  193322,  Aa  and  Abl, 
have  very  broad  and  very  narrow  spectral  lines,  respectively,  and 
indeed  it  is  these  properties  that  can  help  us  distinguish  their 
different  orbital  motions.  We  show  in  Figure  9  CFHT  spectra  of 
the  He  i  7.5876  profile  from  1986  and  2008.  During  this  interval, 
the  broad  component  moved  slightly  redward  as  expected 
for  the  anti-phase  velocity  curve  of  Aa  between  wide  orbit 
phases  0.82  and  0.44  (Figure  8).  We  collected  all  the  available 
red  spectra  that  recorded  He  1 7.5876,  and  we  formed  an  average 
spectrum  for  each  run  in  order  to  increase  the  S /N  of  the  spectra 
at  each  epoch.  We  then  formed  model  spectra  for  each  of  Aa  and 
Abl  from  the  grid  of  Lanz  &  Hubeny  (2003)  using  projected 
rotational  velocities  and  model  parameters  optimized  to  match 


Figure  9.  CFHT  spectroscopy  of  the  He  i  A.5876  line  profile  from  two  epochs. 
The  narrow-lined  component  is  associated  with  Abl,  while  the  broad-lined 
component  corresponds  to  Aa.  The  interstellar  Nai  A.5890  D2  line  appears  near 
Vr  =  +730  km  s_1. 


Table  7 

Radial  Velocity  Measurements  for  Aa 


Date 

(HJD— 2,400,000) 

<t> 

(Wide) 

VT 

(km  s~*) 

O-C 
(km  s"1) 

S/N 

Source 

46608.445 

0.817 

-35.0 

-6.3 

960 

CFHT/ 1986 

49942.734 

0.076 

-39.5 

-40.9 

380 

KPNO/1995 

51467.762 

0.195 

-24.6 

-28.9 

1020 

KPNO/1999 

53246.970 

0.333 

-0.7 

0.3 

240 

OHP/2004 

54675.904 

0.444 

-2.3 

3.8 

750 

CFHT/2008 

55146.400 

0.481 

23.2 

31.1 

220 

NOT/2009 

the  composite  profile  (see  below).  These  model  profiles  were  fit 
to  the  observations  using  a  nonlinear,  least-squares  procedure 
to  derive  the  radial  velocities  of  Aa  and  Abl  at  these  epochs. 
The  derived  Abl  velocities  are  identical  within  errors  to  the 
corresponding  measurements  given  in  Table  5,  and  the  velocities 
for  Aa  are  listed  in  Table  7.  The  errors  associated  with  the 
velocities  of  Aa  are  large,  ±20  km  s-1,  because  this  component 
is  so  broad  and  shallow  and  because  the  shape  of  the  red  wing 
is  sensitive  to  the  details  of  the  removal  of  the  telluric  features 
found  there.  Given  these  larger  errors  and  the  relatively  small 
number  of  measurements,  we  made  a  constrained  fit  of  the 
orbital  radial  velocity  curve  of  Aa  by  setting  all  the  parameters 
from  the  solution  for  Ab  (Table  6)  with  the  exception  of  the 
systemic  velocity  y  and  semiamplitude  K,  and  by  assigning  a 
weight  to  each  observation  proportional  to  the  product  of  the 
spectral  resolving  power  and  the  net  S /N  ratio  in  the  adjoining 
continuum  (Column  5  of  Table  7).  The  fit  (illustrated  in  Figure  8) 
yielded  y  =  — 12  ±  8  km  s-1  and  K  =  21  ±  16  km  s-1,  with  a 
residual  rms  =  21  km  s~  1 .  These  measurements  are  consistent 
with  the  expected  Doppler  shifts  and  masses  for  the  Aa,Ab 
system  (Section  6). 

We  estimated  the  spectroscopic  parameters  for  the  compo¬ 
nents  by  first  reconstructing  their  spectra  using  the  orbital  ve¬ 
locity  curves  and  a  Doppler  tomography  algorithm  (Bagnuolo 
et  al.  1994)  and  then  comparing  the  reconstructions  with  models 
from  the  grids  of  Fanz  &  Hubeny  (2003, 2007).  We  selected  nine 
spectra  from  our  observations  that  recorded  the  blue  portion  of 
the  spectrum  with  a  resolving  power  greater  than  10,000  and  that 
covered  the  extremes  of  motion  in  the  wide  and  close  systems 
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Table  8 

Representative  Stellar  Parameters 


WAVELENGTH  (Angstroms) 


Figure  10.  Example  tomographic  reconstructions  of  the  blue  spectra  of  the 
components  of  HD  193322.  The  solid  lines  show  the  Doppler  tomography 
spectra  while  the  dotted  gray  lines  show  superimposed  model  spectra.  The 
spectra  are  offset  by  steps  of  50%  of  the  continuum  for  clarity. 


(samples  10,  12,  14,  15,  and  16  from  Table  4).  We  began  by  run¬ 
ning  the  tomography  algorithm  for  only  two  components,  Aa  and 
Ab  1 ,  however,  we  found  that  the  subsequent  reconstructed  spec¬ 
trum  for  Aa  had  a  composite  appearance  with  both  broad  and 
narrow  components,  unlike  our  expectation  from  the  He  i  A5876 
profiles  (Figure  9).  We  think  that  this  is  due  to  flux  contamina¬ 
tion  in  our  blue  spectra  from  the  nearby  B  component  (B1.5 
V;  V  sini  ~  100  km  s-1;  see  McKibben  et  al.  1998,  Figure  2; 
Roberts  et  al.  2010,  Figure  1).  Although  component  B  may  be  a 
spectroscopic  binary  with  a  low  semiamplitude  (McKibben  et  al. 
1998),  we  simply  assumed  that  it  was  stationary  and  contributed 
1 1%  of  the  total  flux  (Roberts  et  al.  2010)  in  the  next  iteration  of 
tomographic  reconstruction.  The  power  of  the  tomography  al¬ 
gorithm  to  derive  reliable  and  high-quality  reconstructed  spectra 
increases  with  the  number  of  spectra  and  with  the  orbital  velocity 
range  and  flux  contribution  of  the  components.  Unfortunately, 
in  the  case  of  our  blue  spectra  of  HD  193322,  these  criteria  are 
really  only  met  for  component  Ab  1 .  The  velocity  range  of  Aa, 
for  example,  is  so  small  relative  to  its  characteristic  line  width 
that  the  algorithm  may  incorrectly  assign  line  flux  between  the 
reconstruction  of  Aa  and  the  stationary  component  B.  We  dealt 
with  this  problem  by  starting  the  initial  guess  for  components 
Aa,  Ab2,  and  B  with  model  spectra  rather  than  assuming  a  flat 
continuum  spectrum  (as  done  for  Abl).  Although  the  resulting 
solutions  are  guided  by  our  assumptions,  they  do  at  least  show 
that  the  observed  spectra  are  consistent  with  these  assumptions 
since  otherwise  the  reconstructed  spectra  would  converge  to  an 
appearance  different  from  the  initial  model  guesses. 

We  show  the  results  of  the  full,  four-component,  tomographic 
reconstructions  in  Figure  10.  These  representative  solutions 
were  made  using  the  orbital  solutions  from  Table  6,  adopting  a 
mass  ratio  of  M(Ab2)/M(Abl)  =  0.37  (Section  6)  and  the  flux 
ratios  given  in  Table  8.  These  flux  ratios  were  calculated  from 
the  K'-band  A mciose  and  Amwide  results  (Section  3)  assuming 
that  these  hot  stars  all  contribute  the  same  proportions  to  the 
B  band.  The  spectroscopic  parameters  were  determined  by 
finding  the  Lanz  &  Hubeny  (2003,  2007)  model  that  best 
matched  the  absorption  line  ratios  and  H y  Stark  broadening 
(dependent  on  Teg  and  log  g )  and  with  a  V  sin  i  value  adjusted 
to  fit  the  widths  of  the  absorption  lines  other  than  H y.  The 
results  are  listed  in  Columns  2  and  3  of  Table  8  for  Aa  and 


Parameter 

Aa 

Abl 

Ab2 

B 

F  f  F total  •  ■  • 

0.43 

0.40 

0.06 

0.11 

Teff  (kK)  .  .  . 

33 

32.5 

20 

23 

logg  (cm  s-2) . . . 

4.0 

3.5 

4.0 

4.0 

V  sin  /  (km  s-1) . . . 

350 

40 

200 

100 

M  (  Mq)  . . . 

21 

23 

9 

Mv  (mag)  . . . 

-3.6 

-3.5 

-1.4 

-2.1 

Abl,  respectively,  and  we  estimate  that  the  associated  errors  are 
ATeff  =  ±  1  kK,  Alogg  =  ±0.5,  and  AV  sin;  =  ±40  and 
±10  km  s'1  for  Aa  and  Abl,  respectively.  These  parameters 
suggest  spectral  classifications  of  09  Vnn  and  08.5  III  for  Aa 
and  Abl,  respectively,  based  upon  the  calibration  of  Martins 
et  al.  (2005).  The  “nn”  suffix  for  the  former  classification 
indicates  very  broad  lines.  The  relatively  good  agreement 
between  the  observed  and  model  line  depths  indicates  that  the 
flux  ratios  from  interferometry  (Section  3)  are  fully  consistent 
with  the  derived  strengths  of  the  spectroscopic  features.  The 
parameters  in  Table  8  for  component  B  were  taken  from  the 
work  of  Roberts  et  al.  (2010),  and  the  predicted  model  spectrum 
agrees  well  with  the  narrow,  stationary  spectral  component  from 
the  tomographic  reconstruction.  The  results  for  the  faintest 
component,  Ab2,  are  poorly  constrained  because  this  star 
contributes  such  a  small  fraction  of  the  total  flux,  but  its  spectrum 
suggests  an  early-B,  dwarf  classification.  We  used  the  flux  ratio 
between  Ab2  and  B  and  the  temperature  of  B  and  the  theoretical 
main  sequence  (T^ff,  My)  adopted  by  Roberts  et  al.  (2010)  in 
order  to  estimate  the  effective  temperature  of  Ab2,  assuming  it 
is  a  main-sequence  star.  The  H y  A4340  line  is  the  only  strong 
feature  in  the  reconstructed  spectrum  of  Ab2,  and  the  relative 
weakness  of  the  He  i  AA4387,  447 1  lines  suggests  that  Ab2  may 
also  be  a  rapid  rotator  with  broad  and  shallow  lines.  Note  that 
the  V  sin  i  estimate  for  Ab2  in  Table  8  is  only  approximate  and 
may  be  subject  to  significant  revision. 

6.  DISCUSSION 

One  of  the  primary  goals  of  this  study  was  to  determine 
the  masses  of  the  component  stars.  Since  most  of  our  results 
are  preliminary,  we  cannot  yet  derive  accurate  masses,  but 
the  observational  work  does  demonstrate  the  potential  for  im¬ 
provement  with  further  interferometric  and  spectroscopic  ob¬ 
servations.  The  mass  sums  for  the  wide  and  close  systems 
can  be  determined  from  the  angular  semimajor  axes  and  or¬ 
bital  periods  (Table  2)  and  the  distance  cl  —  741  ±  36  pc 
for  Collinder  419  (Roberts  et  al.  2010).  The  mass  sums 
(Equation  (2))  are  (M(Aa)  +  M( Abl)  +  M(Ab2))/M0  = 
(53.2  ±  11.5)(r//[741  pc])3  forthe  wide  system  and(M(Abl)  + 
M(Ab2))/M0  =  35.6(a/[4.0  mas])3W/[741  pc])3  for  the  close 
system  (where  a  is  the  angular  semimajor  axis  of  Abl,Ab2). 

To  obtain  the  individual  masses,  we  need  to  explore  the 
solution  space  from  the  interferometric  observations  of  the 
close  binary  (Figure  3)  and  from  the  spectroscopic  orbit  of 
Abl  (Table  6).  In  particular,  we  can  use  the  location  of  the 
X2  valley  in  the  (i.  a)  plane  of  Figure  3  to  derive  a  family  of 
solutions  based  solely  upon  a  (since  i  —  i(a)  from  Figure  3). 
We  show  the  derived  individual  masses  as  a  function  of  a  in 
Figure  11.  The  mass  of  Aa  is  set  from  the  difference  of  the 
total  mass  of  Aa,Ab  and  the  mass  of  Abl,Ab2  (from  a  and 
Equation  (2));  the  mass  of  Ab2  is  from  Equation  (3)  (from 
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a  (mas) 


Figure  11.  Diagram  of  the  mass  solutions  for  components  Aa  (solid  line), 
Abl  (dashed  line),  and  Ab2  (dotted  line)  as  a  function  of  a ,  the  angular 
semimajor  axis  of  the  Abl,Ab2  system.  The  shaded  region  surrounding  each 
line  corresponds  to  the  ±la  error  range  in  distance.  The  thick  portion  of 
the  Aa  line  shows  the  section  that  intersects  with  the  d=la  error  range  for 
the  mass  of  Aa  as  determined  from  the  visual  wide  orbit,  assumed  distance, 
and  the  Abl,Ab2  center-of-mass  radial  velocity  curve  (Figure  8).  The  various 
symbols  indicate  the  positions  where  the  mass  ratios  match  those  of  main- 
sequence  stars  with  the  observed  flux  ratios  (squares  for  F{ Ab2) /F( Ab  1 ),  crosses 
for  F(Abl)/F(Aa),  and  diamonds  for  F(Ab2)/F(Aa)).  The  tick  mark  at  the 
bottom  indicates  those  masses  for  which  the  sum  of  the  corresponding  fluxes 
of  main-sequence  stars  attains  a  minimum,  the  situation  most  consistent  with 
the  estimated  total  absolute  magnitude,  and  the  horizontal  line  segment  shows 
the  range  over  which  the  absolute  magnitude  is  within  0.1  mag  of  the  faint 
limit. 

a,  K(Ah  I ),  and  the  i  =  i(a)  relation  in  Figure  3);  and  the 
mass  of  Abl  is  from  the  difference  M(Abl,  Ab2)  —  M(Ab2). 
These  are  all  plotted  in  Figure  1 1  surrounded  by  a  gray  zone 
corresponding  to  the  acceptable  range  in  cluster  distance.  We 
see  that  there  is  a  strict  upper  limit  of  a  <  4.57  mas  required 
to  keep  M( Aa)  >  0.  Furthermore,  we  also  see  that  while  the 
masses  of  Aa  and  Ab  1  cover  a  significant  range,  the  mass  of  Ab2 
changes  little  over  the  range  in  a.  This  is  also  shown  in  Figure  3, 
where  the  location  of  the  y2  valley  is  close  to  a  contour  of 
constant  M(Ab2). 

Another  constraint  on  the  mass  of  Aa  can  be  formed  indepen¬ 
dent  of  the  details  of  the  Ab  1 ,  Ab2  orbit  by  applying  Equation  (3) 
to  the  wide  orbit.  We  take  a,  P,  e,  and  i  from  the  visual  orbit  of 
the  wide  system  (Table  2)  and  combine  these  with  the  orbital 
semiamplitude  K(Ab)  from  spectroscopy  (Table  6)  to  obtain 
M(Aa)/M0  =  (18.7  ±  3.7)(c//[741  pc])2.  The  ±  lcr  region 
from  this  relation  is  plotted  as  the  thick  line  segment  for  M( Aa) 
in  Figure  1 1 ,  and  it  corresponds  to  a  range  in  semimajor  axis  of 
a  —  3.96  ±  0.14  mas. 

It  is  reasonable  to  assume  that  all  three  stars  are  main- 
sequence  objects  given  the  position  of  the  Aa,Ab  system  in 
the  color-magnitude  diagram  (Roberts  et  al.  2010).  The  k'-band 
fluxes  of  massive  main-sequence  stars  scale  with  mass  M  as 
F  oc  M2  30  for  stars  in  this  mass  range  according  to  the  models 
of  Marigo  et  al.  (2008),  so  we  can  use  this  relation  to  predict 
the  flux  ratio  between  any  pair  of  stars  according  to  the  mass 
relations  shown  in  Figure  1 1 .  The  positions  where  the  model 
flux  ratios  match  the  observed  ones  (Table  8)  are  indicated  by 
pairs  of  symbols  in  Figure  1 1 .  These  flux  ratio  relations  indicate 
a  semimajor  axis  range  of  a  =  3.80  ±  0.08  mas. 

A  final  constraint  can  be  set  from  the  overall  fluxes  and 
absolute  magnitude  of  the  combined  system,  ten  Brummelaar 


et  al.  (2000)  estimate  that  the  apparent  E-band  magnitude  of 
Aa,Ab  is  V  =  5.96  ±  0.02  mag,  and  using  the  distance  and 
extinction  for  the  star  from  Roberts  et  al.  (2010),  we  estimate 
that  the  absolute  magnitude  is  My  =  —4.35  ±  0.12  mag. 
The  individual  absolute  magnitudes  of  the  components  are 
given  in  Table  8.  We  can  apply  the  mass-absolute  magnitude 
relation  (M,  My)  for  the  main  sequence  from  the  models 
of  Marigo  et  al.  (2008)  to  obtain  a  prediction  for  the  total 
absolute  magnitude  for  each  of  the  mass  combinations  shown 
in  Figure  1 1 .  The  best  match  occurs  for  the  masses  obtained  at 
a  —  3.81  mas,  approximately  where  Aa  and  Abl  have  equal 
masses.  For  those  masses  the  predicted  absolute  magnitude 
is  My  —  —5.12,  which  is  significantly  brighter  than  the 
estimate  above  from  observations.  The  models  predict  even 
brighter  fluxes  if  either  Aa  or  Abl  is  more  massive,  and  the 
horizontal  line  at  the  bottom  of  Figure  1 1  shows  the  range  in  a 
where  the  combined  magnitude  is  within  0. 1  mag  of  the  faint 
limit. 

The  average  estimate  for  the  semimajor  axis  from  the  above 
three  constraints  is  a  =  3.85  ±  0.09  mas,  and  we  adopt 
the  associated  mass  solution  as  best  representing  the  current 
observational  data.  The  masses  and  other  properties  summarized 
in  Table  8  are  generally  in  agreement  with  expectations  for 
hot,  main-sequence  stars  (Martins  et  al.  2005).  However,  a 
number  of  significant  discrepancies  remain  that  deserve  further 
investigation.  The  mass  of  component  Abl  is  similar  to  that 
expected  for  an  08.5  III  star  («24M0;  Martins  et  al.  2005), 
but  the  star’s  absolute  magnitude  is  about  1.8  mag  fainter  than 
typical  for  such  stars.  This  discrepancy  hints  that  Abl  may  be  a 
dwarf  rather  than  a  giant  star.  The  overall  faintness  of  the  system 
compared  with  expectations  for  the  stars’  masses  may  indicate 
that  the  distance  estimate  needs  to  be  revised  downward  (leading 
to  lower  masses)  and/or  the  extinction  estimate  revised  upward. 
Some  confusion  also  remains  about  whether  Aa  or  Abl,Ab2 
is  brighter.  As  we  noted  in  Section  3,  the  visibility  analysis 
indicates  that  Amwide  =  m(Aa)  —  m(Ab)  =  0.086  ±  0.012  mag 
(Ab  brighter  than  Aa),  which  agrees  within  errors  with  high 
angular  resolution  measurements  by  Malz  Apellaniz  (2010) 
using  the  AstraLux  camera,  m(Aa)— m(Ab)  =  0.04  ±  0. 19  mag. 
Although  these  results  indicate  that  Ab  is  somewhat  brighter 
than  Aa,  we  refrain  from  re-designating  the  identities  of  the 
components  to  avoid  confusion  with  published  results. 

We  expect  that  these  lingering  problems  will  be  resolved 
with  future  interferometric  observations  that  will  better  sample 
the  orbital  phases  of  the  close  pair  near  the  nodal  crossings 
and  will  lead  to  improved  constraints  on  the  angular  semimajor 
axis  a.  In  addition,  we  clearly  need  to  continue  the  long-term 
high-resolution  work  on  the  wide  orbit  to  cover  the  missing 
orbital  phases  (Figure  2).  We  plan  to  obtain  these  measure¬ 
ments  through  continuing  observations  of  this  system  using  the 
CHARA  Array  interferometer.  Additional  high  S/N  and  high- 
resolution  spectroscopy  holds  the  promise  of  delivering  better 
orbital  constraints  on  Aa  and  Ab2  that  would  then  allow  us  to  es¬ 
timate  the  masses  without  relying  on  the  distance  of  the  cluster. 
Indeed,  reliable  orbital  elements  would  render  it  possible  to  set 
an  independent  estimate  of  the  cluster  distance.  We  anticipate 
expanding  the  spectroscopic  coverage  over  the  next  decade. 

We  close  with  some  speculative  remarks  about  the  angular 
momentum  distribution  of  the  stars  of  HD  193322.  It  is  re¬ 
markable  that  this  system  contains  both  a  very  rapidly  rotating 
star  (Aa)  and  a  very  slowly  rotating  star  (Abl).  It  is  possi¬ 
ble  that  Abl  has  a  rotational  axis  with  a  low  inclination,  so 
that  its  equatorial  rotational  velocity  is  close  to  typical  values. 
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However,  the  very  large  line  broadening  of  Aa  places  it  among 
the  most  rapidly  rotating  O-type  stars  known  (Penny  1996).  It 
is  possible  that  the  angular  momentum  of  the  natal  cloud  led 
directly  to  rapid  rotation  in  the  case  of  Aa  and  to  the  forma¬ 
tion  of  a  binary  in  the  case  of  Ab.  Alternatively,  there  may 
have  been  some  very  close  gravitational  encounters  in  the  early 
life  of  the  system.  In  some  circumstances,  a  close  encounter 
between  a  binary  and  a  third  interloper  can  lead  to  a  merger 
of  two  of  the  components  and  ejection  of  the  third  (Gaburov 
et  al.  2010).  It  is  possible  that  the  rapid  rotator  Aa  is  such  a 
merger  product  and  that  the  runaway  star  68  Cygni  was  the  ob¬ 
ject  ejected  from  the  system  (Schilbach  &  Roser  2008).  If  so, 
then  the  orbital  and  spin  properties  of  the  stars  of  HD  193322 
offer  key  evidence  about  the  early  dynamical  processes  in  this 
cluster. 

The  CHARA  Array,  operated  by  Georgia  State  University, 
was  built  with  funding  provided  by  the  National  Science  Foun¬ 
dation,  Georgia  State  University,  the  W.  M.  Keck  Foundation, 
and  the  David  and  Lucile  Packard  Foundation.  This  material  is 
based  upon  work  supported  by  the  National  Science  Founda¬ 
tion  under  grants  AST-0606861,  AST-0606958,  AST-0908253, 
and  AST- 1009080.  Institutional  support  has  been  provided  from 
the  GSU  College  of  Arts  and  Sciences  and  from  the  Research 
Program  Enhancement  fund  of  the  Board  of  Regents  of  the 
University  System  of  Georgia,  administered  through  the  GSU 
Office  of  the  Vice  President  for  Research.  B.D.M.  and  W.I.H. 
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